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Abstract. The magnetic pyrochlore oxide materials of general chemical formula 
Pt2Ti 2 07 and R^S^Oy (R = rare earth) display a host of interesting physical 
behaviours depending on the flavour of rare earth ion. These properties depend 
on the value of the total magnetic moment, the crystal field interactions at each 
rare earth site and the complex interplay between magnetic exchange and long- 
range dipole-dipole interactions. This work focuses on the low temperature physics 
of the dipolar isotropic frustrated antiferromagnetic pyrochlore materials. Candidate 
magnetic ground states are numerically determined at zero temperature and the role of 
quantum spin fluctuations around these states are studied using a Holstein-Primakoff 
spin wave expansion to order 1/S. The results indicate the strong stability of the 
proposed classical ground states against quantum fluctuations. The inclusion of long 
range dipole interactions causes a restoration of symmetry and a suppression of the 
observed anisotropy gap leading to an increase in quantum fluctuations in the ground 
state when compared to a model with truncated dipole interactions. The system retains 
most of its classical character and there is little deviation from the fully ordered moment 
at zero temperature. 
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1. Introduction 

Magnetic frustration arises when magnetic moments (spins) are unable to minimize 
their classical ground state energy by minimizing the two body magnetic exchange 
interactions pair by pair. The simplest example of magnetic frustration is for 
antiferromagnetically coupled spins on a triangular lattice where the triangular 
coordination prevents the spins from pointing mutually antiparallel to each other and 
where, consequentially, the spins adopt a non-collinear structure to minimize their 
energy Due to the canted structure, the internal mean-field is much smaller in 
frustrated spin systems than in conventional collinear magnets. As a result, frustrated 
magnetic systems have an enhanced propensity for large zero temperature quantum spin 
fluctuations allowing for the possibility of exotic quantum ground states [TJ |2] . 

The pyrochlore lattice of corner- sharing tetrahedra where spins interact via nearest 
neighbour antiferromagnetic exchange interactions has attracted much interest over the 
past fifteen years [3- In particular, the R,2Ti 2 07 and R2Sn 2 07 magnetic insulators 
(where R is a trivalent rare-earth ion, R=Tb, Ho, Dy, Gd, Yb and Er, sitting on a regular 
pyrochlore lattice) have been found to display a multitude of magnetic phenomena. For 
example, the combination of the strong Ising anisotropy in Dy2Ti20"7 jHEj, Ho 2 Ti207 
[HI Cj) an d Ho 2 Sn 2 7 [H] gives rise to "spin ice" behavior where the system exhibits 
extensive magnetic entropy and maps onto the problem of proton disorder in common 
ice water [5J (TUJ . In Yb 2 Ti 2 07, a sharp first order transition is observed in specific 
heat measurements, but with no obvious sign of long range order in neutron scattering 
The Tb2Ti 2 07 antiferromagnet is particularly interesting since it appears to be 
an excellent candidate for collective paramagnetism ^2] or "spin liquid" behavior, as 
it fails to develop long range order down to 50 mK despite a Curie- Weiss temperature 
of CW = -20 K [13 d EH- Finally, both Er 2 Ti 2 7 PU and Gd 2 Ti 2 7 [13 dl d] 
display a transition to long-range ordered states below 1 K. 

In the Ho 2 Ti 2 07, Dy 2 Ti 2 7 , Tb 2 Ti 2 7 and Gd 2 Ti 2 7 materials, the rare-earth 
ion carries a large magnetic moment between five and ten Bohr magnetons and the 
dipolar energy scale at nearest neighbour makes a sizeable contribution to the total 
magnetic interaction. Indeed, in Dy 2 Ti 2 07 and Ho2Ti 2 07, it has been shown that it 
is the dipolar interactions that cause the spin ice phenomenology jSlE]- In Tb2Ti 2 07, 
the dipolar interactions play a very important role in frustrating the nearest neighbour 
exchange interactions and are a significant player in the formation of a possible spin 
liquid state In Gd2Ti 2 07, the 4/ orbital of Gd 3+ is half filled, and crystal-field 

single ion anisotropy effects are negligible. This system, with #cw ~ -9.8 K, and spin 
quantum number S = 7/2 therefore constitutes a useful starting point as a real material 
which could be described by a nearest neighbour Heisenberg pyrochlore model ffi I20j. 
However, in this system the large dipolar interactions are important since they have 
been predicted to select a unique long-range ordered classical ground state [THl 1211 

Rare-earth pyrochlores do not hold a monopoly on unusual phenomena among the 
insulating rare-earth magnetic materials. Possibly one of the most interesting cases is 
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that of the GdsGasO^ three dimensional garnet of corner- sharing triangles [22j. In this 
material, experimental evidence of a magnetic field induced spin liquid has been found 
[2*2*] . as well as reentrance from a long range ordered magnetic state to a spin liquid 
upon cooling [231 121! • 

The discussion in the previous two paragraphs brings up the following question. In 
conventional exchange-coupled Heisenberg antiferromagnets the quantum fluctuations, 
AS, are governed by an integral over the Brillouin zone of l/u(k) where u(k) are 
the fe-dependent excitation frequencies. Weakly dispersive modes give rise to larger 
fluctuations and (anisotropy) gapped modes contribute smaller quantum fluctuations 
|25j . In the nearest neighbour Heisenberg pyrochlore, we have two zero energy modes 
throughout the zone and any conventional perturbative spin wave calculation is therefore 
singular. In systems such as the rare-earth pyrochlores above, we expect both single ion 
anisotropy and dipolar interactions to lift the classical degeneracy and make the spin 
excitations dispersive. The leading effect of those perturbations will be to push the two 
zero modes throughout the zone to finite frequency. This is akin to the finite energy 
fc-independent resonance observed in inelastic neutron scattering in the ZnCr 2 04 spinel 
[2E! . In other words, in the isotropic Heisenberg antiferromagnet with conventional long 
range Neel order, it is the dispersion of the acoustic magnons {uj(k) ~ |fe|) that controls 
the leading 1/S quantum fluctuations. The question therefore concerns the magnitude 
of quantum fluctuations arising from gapped excitations with negligible dispersion 
throughout the Brillouin zone in a highly frustrated system of the type considered here. 
This would seem particularly relevant to Gd2Ti2C>7 where the experimentally observed 
ground state seems to be at odds with classical predictions 121] • 

The outline of the paper is as follows. In Section 2 we discuss the general framework 
for calculating the quantum spin fluctuations in a non-Bravais lattice in the presence 
of combined exchange and long range dipolar interactions. In Section 3, we apply the 
method discussed in Section 2 to the specific case of the Gd2Ti 2 07 pyrochlore material. 
Section 4 contains a brief conclusion of our work. 

2. 1/S Formalism 

In this section we introduce the Hamiltonian believed to describe our antiferromagnetic 
pyrochlore system and map out its low temperature excitations using linear spin wave 
theory. From the diagonalized form of the spin wave Hamiltonian we obtain expressions 
for the zero point energy shift and reduction in the classical staggered moment due 
to spin fluctuations. Utilizing the partition function of the non-interacting Bose gas, 
the relevant thermodynamic quantities are given in terms of the calculated spin wave 
dispersion energies. 
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2.1. Model Hamiltonian 

A general two-body spin interaction Hamiltonian including both p th nearest neighbour 
anisotropic exchange and long range dipole-dipole interactions can be written as 



H = -\ E E E J S i R 5) s k a (Rns i ,(Rn 

's a (R») ■ S P {R>) (W • fl£5) ^(Jf ) • 



2 



+ ^EE 



using flg = (JF + r/J ) - + r a ) , J*j, (J^) = En=i Jnf^n^-RZ,, is the set 

of n th nearest neighbour connection vectors) and D^d = A i o(5'A i B) 2 /4vr where the spins 
are full 0(3) operators satisfying S Q (-R M ) • S^R^) = S(S + 1). A factor of | has been 
included to avoid double counting. The anisotropic exchange tensor J® is the relative 
strength of the magnetic interaction for spin components k and I between n th nearest 
neighbours. The factor f^n is equal to 1 — 5 a ,p if the coupled spins are on different 
sublattices, or 5 at p if they are on the same sublattice. The symbols \x and v index all 
primitive Bravais lattice vectors, a and f3 run over all basis vectors and k and / over the 
spin components x,y and z. We wish to consider the stability of some classical ground 
state by investigating the role of quantum fluctuations in reducing the fully polarized 
classical spin value S = (0,0, S). This can be accomplished by changing the axis of 
quantization from the global ^-direction (which is an arbitrary choice of the theory) 
to a local quantization axis which points in the direction of each spin. Let S a (R^) 
be the locally quantized spin which is related to the real spin operator defined in the 
crystallographic frame via the rotation 

S a (R! t ) = R(e a ,<l> a )S a (R! t ). (2) 

This transformation can be implemented in an alternate way by defining an orthogonal 
triad of unit vectors n i a for each sublattice a by 

n i>Q = R(9 a ,(f>a)ei (3) 

where {e x ,e y ,e z } are the usual Cartesian unit vectors and clearly n za points in the 
direction of the spin S^R 11 ). In this form, we can re-express Equation (J2j) as 

S a (R^) = J2s i a (Rnn i , a (4) 

i 

and substituting into our original Hamiltonian Equation (JIJ we find, 

H = ~\ E E E z% ( R Z) WW (5) 

a,/3 i,j /J.,u 



where 



k.l 



(6) 
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is the two spin coupling matrix. Making use of the periodicity of the lattice we Fourier 
transform the spin operators via 



fc 



where iV is the number of Bravais lattice points, leading to 



2 

o,/3 ij n,v fc,fc' 



x exp[ifc • (RP + r a )] expfifc' ■ (RY + r^)] (8) 

and using the definition of the Dirac delta function 8 qq i = -k ^„exp[i(qr — q') ■ R v ], we 
can write this as 

n = ~\ E E E ( fc ) ^(- fc ) ( 9 ) 



2 

a,/3 i,j fc 



where the Fourier transform of the spin-spin interaction matrix has been defined by 

3% ( k ) = J2 4* ( R Z) exphifc " (10) 



ft 

The Holstein-Primakoff transformation |2Zj defined by 



S+iRP) = A/2S - cl(R»)c a (R») c a {R») (11a) 



S~(R») = cUR") V 2S - cL(R»)c a (R») (lib) 

S z a (R») = S-cl(Rnc a (Rn (11c) 

allows us to express our anticommuting fermionic spin operators in terms of commuting 
bosonic spin deviation annihilation and creation operators c a (R^) and c^(J2 M ) where 
[c a (J?^), cL(R u )] = SafiSfrv. We wish to deal with a theory whose transverse terms 
are only linear in the boson operators, leading to a quadratic boson-mode spin wave 
Hamiltonian. In other words, a harmonic theory where there exists only small deviations 
away from some stable minimum energy spin configuration. Identifying that the boson 
operators can be written in terms of their Fourier components 

c a (R») = -L c -( k ) expRfc • (R^ 1 + r«)] (12a) 
fc 

4 = 4= E c -( fe ) ex pt ifc • + r «)l> ( 126 ) 

V ^ fc 

we can write the linearized Holstein-Primakoff transformation in reciprocal space as 

S X M = ^ [cl(k) + c a (-k)] + O f (13a) 

~ S i( k ) = [cj, (fc) - c a (-fe)] + O (j-) (135) 
S z a (k) = VNS6 kfi exp[-ik-r a }--^=J2 c «( k ') c <*( k '- k )> ( 13c ) 



N fc' 
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where the Fourier transformed boson operators satisfy [c a (k), ct(fc')] = b~ a> pb~k,k'- 

Performing this transformation on Equation © and keeping only terms up to 
second order in the boson mode amplitudes, we obtain a form for the Hamiltonian 
to order 1/S which can be broken into three parts that are zeroth, first and second 
order in the boson operators. 

n = n {0) + n {1) + n {2) (14) 

with 

H (o) = _ Ijvsp + 1)5^(0) (15a) 

= "Sa/^E [ F ^( )4(0) + F^(0)c a (0)] (156) 

a,/3 
a,/3 k 

+ B* ap (-k)c a (-k)cp(k) + A; /3 (-fe)c a (-fc)4(-fe)] (15c) 

and 

F Q/3 (0) =J^(0)+iJ"5(0) (16a) 

A^(fc) = \ {J X J (*) + J$ (k) - i [JS (k) - J$ (k)} } - J%Waj> (166) 

7 

B^(fc) = \ (k) ~ J$ (k) + i [JS (k) + J$ (*)] } • (16c) 

Here we note the manifestation of the equilibrium condition, that for a particular 
classical ground state to be stable, the effective magnetic field acting on each spin is 
parallel to that spin. This implies that there are no terms linear in the boson operators in 
Equation JUJ), or J2p M ) = showing that the J% (0) , J% (0) , J$ (0) and J% (0) 
blocks of the interaction matrix are only responsible for determining the equilibrium 
configuration and do not play any dynamical role |28j . 

2.2. Diagonalization 

If a stable ground state spin-configuration exists, then the quadratic spin wave 
Hamiltonian in Equation ()15cj) can be rewritten in terms of a general quadratic form 
by defining a vector which contains the boson creation and destruction operators 

Cl (k)...c s (k), c \(-k)...cl(-k)} T (17) 



X 



where s is the number of atoms in the basis. In matrix notation, this leads to the 
Hamiltonian 

n = n {0) + ^x f Hx (is) 
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where H is the 2s x 2s matrix defined by 



H 



is 

2 



A(fe) B(fc) 
B*(-fe) A*(-fe) 



(19) 



and A(fe) and B(fc) are the full matrix representations of Equation (jl6fc|) and (|16cj) . 
We now introduce the transformation X = QY where Y is a column matrix of new 
boson operators a a (k) and a) a (k). It is clear that if we intend to diagonalize our spin 
wave Hamiltonian, the transformation law must satisfy Q^HQ = A where A is the 
diagonal matrix of eigenenergies. The matrix Q constructed from the eigenvectors of 
the transformation must be normalized with respect to a metric which preserves proper 
boson commutation rules for our new operators. This requirement leads to the possibility 
of a non-unitary Q, and in general forces us to perform a Bogoliubov transformation to 
diagonalize the spin wave Hamiltonian The details of the numerical construction 
of this transformation are included in |Appendix A| and the result is given by 



W = H<»> + yjyj £ „(fc) 



EE* 



(fe) [o) a {k)a a {k) + al(—k)a a (—k 



(20) 



2.3. Zero Point Effects 

In a quantum mechanical antiferro magnet, there exists spin fluctuations at zero 
temperature with energy e(k) which raise the classical ground state energy and reduce 
the staggered magnetic moment per spin from its classical value of S. Examining 
Equation (f2T)|) and noting the fact that at zero temperature, (a£(fe)a a (fc)) = 0, the 
contribution to the ground state energy from quantum spin fluctuations to order 1/S 
will be given by 

A ^ (0) = EE^ fc )- ( 21 ) 

fc a 

Returning to the original Holstein-Primakoff bosons, the reduction in the classical 
spin polarization will be equivalent to the average number of activated magnons localized 
at each site, 

AS = ^EE<^) c ^)> ( 22 ) 

fi a 

where (• • •) indicates the expectation value with respect to the classical ground state. 
Upon moving to Fourier space, this becomes an average over the first Brillouin zone 
given by 

AS = ^ E E + 4(-fcM-*)> = 2^ E < xtx > - \ (23) 

fc a fc 

Recall that we have introduced new boson operators a a (k) which diagonalize the spin 
wave Hamiltonian through the canonical Bogoliubov transformation X = QY, where 
it is convenient to note that the transformation can be written in component form as 
X a = J2i Qa.iYj and X Q+S = Q a+s ^Yi with the convention that Greek indices run 
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between 1 . . . s and Latin indices between 1 ... 2s. The inner product, X^X can then be 
expanded as 

XtX = EE [^ a ,^ a {k)ap{k) + Ql +Sti Q^ +s al(-k)a a (-k) 

i a,/3 

+ Q^Q^ +S 4 (fe)4(-fe) + QU.iQW"*)^*) + QUA,/3 +s <W] (24) 

and at zero temperature there are no thermally activated spin deviations to destroy with 
a a (k), so the expectation value of all quadratic boson terms tend to zero. Substituting 
Equation into and neglecting these terms we find the staggered moment at 
zero temperature due to quantum spin fluctuations to be 



\ k a / 



AS = - — VV^Q] -1 • (25) 



At this point we may reflect upon the fact that beginning with Equation (fT%|) the 
form of the fully diagonalized spin wave Hamiltonian to 0(1/S) in Equation pUj) could 
have been inferred with little or no calculation. Physically, to order 1/S, the Hamiltonian 
ought to describe a system of non-interacting bosons with energy dispersion equal to 
Huj a (k) where oJ a {k) are the classical spin wave frequencies obtained by linearizing the 
classical equations of motion for interacting magnetic dipoles. Self contained approaches 
exist which begin from such a Hamiltonian and allow for the calculation of the staggered 
magnetic moment at zero temperature without any knowledge of the eigenvectors of 
the diagonalizing transformation embedded within the Q matrix discussed above. A 
staggered anisotropy field h a is introduced which points along the local quantization 
(classical ground state) direction and the change in the ground state energy dE/dh a is 
calculated giving the total magnetic moment on each sublattice at T = K [30J. In the 
discussion above, we have chosen not to follow this physically motivated approach and 
have instead utilized the full spectrum of our transformation matrix Q to compute the 
desired zero-point quantum effects. 

2.4- Thermodynamic Relations 

We have decomposed our theory of spins on a lattice interacting via frustrated exchange 
and dipole-dipole interactions to that of non-interacting Bose gas (Equation (j20J))- 
Beginning from the well known partition function, 

Z = Vevl-m = exph/TC] g g l _ exp[-fe, W ] (26) 

where E = H {0) + AH (0) , and is the inverse temperature we have used the standard 
thermodynamic relations to derive equations for the free energy (JF), average internal 
energy (U), specific heat at constant volume (C v ) and the entropy (S) as a function of 
temperature. These expressions will be valid only for low temperatures where T e a (k) 
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and are given by 

T = E + /T 1 E ln - ex P[/M fc )D U = E + E £a(fe)n B (e a (fc)) 

k a k a (27) 

C - = /3 2 ^^[^(fc)n B (e a (fc))] 2 exp[^ a (fc)] S = P{U-T) 

k a 

where n B (e a (k)) = (exp[f3e a (k)] — is the Bose factor [31. In order to calculate 
the magnetic order parameter as a function of temperature (m*), we must return to the 
general form for the reduction of the sublattice magnetization in Equation ()24j) . For 
suitably low temperatures, the expectation value of all non-number conserving quadratic 
boson terms will be negligible and we find 

m* = S - AS(T = 0) - ^EE [tfQL^Mfc)). (28) 

k a 

3. The Heisenberg-Like Dipolar Pyrochlore 

Having developed the formalism of linear spin wave theory for a multi-sublattice non- 
Bravais magnetic system we now map out the low temperature excitations of spins 
interacting via nearest neighbour Heisenberg exchange and long range dipole-dipole 
interactions on the pyrochlore lattice. 



3.1. Classical Ground States 

It is known that when isotropic nearest neighbour exchange interactions alone are 
considered on the pyrochlore lattice, the system exhibits no global symmetry breaking 
and has soft modes at long wavelengths p3 1201- The inclusion of a perturbation such 
as long-range dipole-dipole interactions breaks the continuous 0(3) symmetry [T%t I2T] 
and leads to a gap in the excitation spectrum. The non-trivial coupling of the isotropic 
spin symmetry and translational lattice symmetry for a particular set of interaction 
parameters controls how these optical modes disperse. The details of this dispersion 
directly affect the stiffness of a given spin configuration and quantitatively determines 
the level of proliferation of quantum fluctuations in the ground state. 

Palmer and Chalker have investigated the classical Heisenberg model for spins 
residing on the pyrochlore lattice with long-range dipolar interactions |21j. They find 
that the dipolar perturbation fixes the two continuous internal degrees of freedom and 
leads to the energetic selection of the six degenerate q = (Palmer-Chalker, PC) ground 
states shown in Figure ^ where the total magnetic moment on each tetrahedron is 
zero. We have confirmed the six degenerate PC classical ground states numerically in 
real space using a field averaging approach by exploring the energy landscape at zero 
temperature as shown in Table ^ The inclusion of long range dipole interactions does 
indeed break the global rotational symmetry and selects a crystallographic plane in 
which the spins are confined while still preserving a zero moment structure on each 
tetrahedron. The discrete degeneracy of the PC ground states allows us to perform our 
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Figure 1. Six degenerate Palmer-Chalker (PC) ground states (reverse spins for other 
three) for spins on a single tetrahedron. 

Table 1. Vectors which define the six classical PC ground state spin configurations 
on the pyrochlore lattice. 





State 


PCxz 


State 


PCyz 


State 


Si = 


1/V2(1,T,0) 


Si = 


1/^2(1,0,1) 


Si = 


1/^2(0,1,1) 


s 2 = 


1/^2(1,1,0) 


S-2 = 


1/^2(1,0,1) 


s 2 = 


1/^2(0,1,1) 


S3 = 


1/^2(1,1,0) 


s 3 = 


1/^2(1,0,1) 


S3 = 


1/^2(0,1,1) 


s 4 = 


1/^2(1,1,0) 


s 4 = 


1/^2(1,0,1) 


s 4 = 


1/V2(0,1,1) 



spin wave calculation on only one state, and the remainder of the results presented 
were calculated for the PC xy state, but are identical (with the proper symmetry 
considerations) for the other five states. 

3.2. Spin Wave Dispersion 

The functional dependence of the normal modes of an oscillating system on the 
wavevector can be derived from classical mechanics. Treating each spin as a classical 
rotor, the rate of change of the angular momentum in the local spin frame, h~dS a (R^) / dt 
is equal to the torque, N a (R^) = pL a {R^) x B a (R^) acting on 5 a (i? M ), where B a (W l ) 
is the total effective magnetic field experienced by a spin on site a in the local frame. 
Here we must take care to ensure that the vector product is carried out in a way 
that is independent of the local coordinate system. Following this procedure we use a 
plane- wave ansatz of the form S^R^) ~ exp [ik ■ (i? M + r a ) — iut] to solve the resulting 
system of linear equations and obtain the excitation spectrum. The spin wave dispersion 
spectrum calculated using Holstein-Primakoff bosons above can then be compared with 
the classical modes of excitation as an essential consistency check. 

Armed with the knowledge of a set of candidate classical ground state spin 
configurations for the isotropic dipolar pyrochlore (which satisfy Ti^ = 0) we may 
now attempt to apply this method to an actual material, the frustrated pyrochlore 
antiferromagnet gadolinium titanate. Gd2Ti 2 0y is a spin only (S = 7/2) material with 
an isotropic magnetic exchange parameter between nearest neighbours which can be 
estimated from the value of the Curie- Weiss constant (#cw = —9.8 K) using the relation 
J\ = 36 l cw /[zS(S + 1)] ~ —0.30476 K where z is the number of nearest neighbours. The 
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dipolar energy scale is set at D M /(ay/2/A) 3 ~ 0.05338 K where a = 10.184 A is the edge 
size of the cubic unit cell and ITZ 1 ] = ay/2 /A is the nearest neighbour distance leading 
to (AalttWlJil ~ 1/5 P2j. 

In order to map out the low temperature magnetic excitations of this material 
we are faced with the non-trivial task of calculating the Fourier transform of the 
12 x 12 interaction matrix J% (Equation ()15cj) ) which contains the exchange 

contribution out to a finite number of nearest neighbours, and the infinite range dipole- 
dipole interactions. This is accomplished via the Ewald summation method [32] and the 



final result for up to third nearest neighbour isotropic exchange 
3% (*) = ^ (*) - D dd V% (fc) 



is given by 



where 

'a/3 



El 3 a (k) 



and 



U al3 



(fc) 



(29) 



2J 1 (n i>a • n jt p) cos[fc • (r a - r /3 )](l - 5 a>p ) 
+ 2J 2 (n iiQ ■ n ji/3 ) ^ cos l k ' ( r « +r/3- 2r 7 )](l 

+ 2J 3 ( )y^cos[2fc ■ (r CT - r p )]S, 



(30) 



cr<p 



4A 3 



3^ 



n^)<W + ^ £ /Cg, (fc - G) e"l fc - G lV^ 2 e xp[-iG ■ r Q/3 ] 



+ [Ml* (flg) - W2« (Jig)] exp[-ifc ■ Jig 



(31) 



with 



^3, (fc 



G) 



[n i>a • (fc - G)] [n hP • (fc - G)] 



in,- 



™i,/3) 



fc - G| 
1 

I n^ v 1 2 
I ^ | 



2A 



7T 



-A 2 R% 

I a P I 



+ 



erfc(A|i?T 



'a/31 



X 



2 A 



2A 2 + 



1 1?^" I 

I ""a/? I 



-a 2 Lre 



a/3 



a OTfc(A|flg| 



I F?^ I 



(32a) 
(326) 



(32 c) 



In these expressions A = \/tt/v is the Ewald convergence parameter, erfc(- • ■) is the 
complimentary error function, v is the volume of the primitive unit cell, G is a FCC 
reciprocal lattice vector and the primed summation indicates that the real-space sum is 
over all FCC Bravais lattice vectors except R^L = 0. 

If we completely ignore the effect of dipole interactions we expect the excitation 
spectrum of a system with massless Goldstone bosons as seen in Figure El (A). There are 
two unfixed internal degrees of freedom associated with the spin isotropy and one time 
reversal degree of freedom. The ground state spin configuration can fluctuate with no 
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Figure 2. Spin-wave energies, e a (k) in kelvin, for (A) the pure Heisenberg pyrochlore 
{Ddd = 0) with Ji = —0.30476 K where both branches are doubly degenerate and (B) 
excitations from one of six discretely degenerate classical ground states for Gd 2 Ti 2 07 
with nearest neighbour exchange and both nearest neighbour and long range dipole 
interactions using the Ewald method. 

cost in energy between states in its continuously degenerate manifold while maintaining 
a macroscopic zero-moment structure. 

When the effects of dipole-dipole interactions are included, the internal degrees of 
freedom are fixed, and the spins are able to minimize their energy by ordering in a non- 
collinear, planar configuration which satisfies a pair- wise antiferromagnet constraint. 
This can be done in a number of ways, provided that each spin is parallel to an opposite 
edge of the tetrahedron (refer to Figure [TJ and we move from a continuously degenerate 
manifold to one with the six Palmer- Chalker discrete ground states. 

Beginning from the PC xy state with isotropic nearest neighbour exchange {J 2 = 
J 3 = 0) the spin wave dispersion spectrum for the parameter set of Gd 2 Ti 2 7 is shown 
in Figure El (B). We observe four optical modes indicating that all distortions from 
the chosen ground state are associated with a positive energy cost. As predicted, the 
introduction of anisotropy has broken a global symmetry and lifted the spectrum to finite 
frequency. The remnant soft modes from the pure Heisenberg model (Figure El (A)) have 
remained moderately degenerate and relatively dispersionless, while the acoustic modes 
have split and are dispersing to the point where one mode has changed curvature along 
r — ► X (for the PC xy state). Once the system is in one of its stable classical ground 
states, it is unable to move between equivalent energy states through a low energy 
long-wavelength acoustic mode. This implies that we were justified in making our semi- 
classical approximations as it is likely the Gd2Ti207 may display very little quantum 
mechanical behaviour to asymptotically low temperatures. 

The inclusion of truly long range dipole interactions using the Ewald technique, as 
opposed to truncating at nearest neighbour, leads to a restoration of some previously 
buried lattice symmetry. In Figure El (B) we observe promotion of quasi-degeneracy in 
the excitation spectrum between W and K, and a slight suppression of the anisotropy 
gap throughout the Brillouin zone due to the global dipolar field being partially screened. 
It is interesting to note that when interaction terms beyond nearest neighbour are 
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neglected, quantum fluctuations are inhibited in the ground state as the local unscreened 
interaction leads to an increased anisotropy gap at the zone center; e(fe) using the Ewald 
method is lower than s(k) using only nearest neighbour interactions at V. 

3.3. Quantum Spin Fluctuations 

To form a quantitative picture of the role of spin fluctuations in the ground state we 
must calculate both the shift of the classical ground state energy and the reduction of the 
fully polarized S = 7/2 moment in Gd2Ti 2 07 at zero temperature. This is accomplished 
using the properly normalized eigenvectors of the numerically constructed Bogoliubov 
transformation in conjunction with Equations ()15aJ) and (}2~Tj) for the energy of the system 
and Equation (j25)) for the reduced moment. Performing averages over the first Brillouin 
zone of the FCC lattice, the main quantitative numerical results of this study, accurate 
to order 1/S are 



where E = — 1/2NS 2 J2 a a <?ap (0) * s the classical ground state energy, NN = 
Nearest Neighbour Dipole and LR = Long Range Dipole. In accordance with our 
naive expectation from an examination of the spin wave dispersion spectrum, we see 
a softening of the system and enhanced spin fluctuations in the ground state when 
long range dipole interactions are treated properly. The reduction of the fully polarized 
moment is almost twenty percent larger when the Ewald summation method is employed. 
These comparisons are overshadowed by the reality that the actual numerical values of 
the quantum corrections are negligible when the zero point shifts above are divided 
by the large S = 7/2 moment of the Gd 3+ ion. The ground state energy is shifted 
by less than ten percent and the fully ordered moment is reduced by half of that. 
However, the S-dependence of the numerical values for the quantum corrections in 
Equation (j33a|) and (J33fej) has been explicitly removed, and these numbers are set by J 
and (;Uo/47r)((yf/iB) 2 / \TV~^ for a given interacting system. It is apparent that one can 
not a priori rule out the possibility of sizeable quantum fluctuations (~5%forS = 7/2 
here) in the ground state of the frustrated pyrochlore materials. Yb2Ti207 falls into 
the category of an XY S = 1/2 system [TT] and using Equation ()33aj) and ()33fc|) above 
we might expect quantum fluctuations on the order of 20%. Therefore, it is crucial to 
properly identify and include the variety of complex competing interactions which act 
as finite, non-zero perturbations to the isotropic Heisenberg exchange parameter J. 

3.4- Thermodynamic Properties 

We have calculated the specific heat and staggered magnetization at finite temperatures 
and the results are shown in Figure El At low temperatures, all four spin wave modes are 
gapped and we observe an exponentially vanishing specific heat as T — > K unlike the 
T 2 laws observed for both Gd 2 Ti 2 07 (HI and Gd 2 Sn 2 7 j^Zj. It appears that the energy 



E NN = E (l + 0.1442/S) 
E LR = E (l + 0.3526/S) 



m* NN = S(l - 0.1007/S) 
m* LR = S(l - 0.1174/S) 



(33 a) 
(335) 
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T [K] T [K] 

Figure 3. The specific heat at constant volume scaled by the temperature (A) and 
the staggered magnetic order parameter (B) as a function of temperature calculated 
using linear spin wave theory for gadolinium titanate. 

scale set by the anisotropy gap corresponds to the peak position near 0.4 K in panel (A) 
and there exists some qualitative similarities with the experimentally measured value of 
C v /T [TS] where two peaks are observed at approximately 0.7 K and 1.0 K. However, 
we have confined the system to a particular classical ground state and cannot probe the 
existence of multiple phase transitions in this study [THl EH EH] • 

The temperature dependence of the staggered magnetization shown in Figure El 
panel (B) is also of interest and we find the behaviour, 1 — exp[— A/T]T 3 , with A ~ gap, 
as expected for an antiferromagnet at low temperature [28] . 

4. Conclusions 

We began by examining the energy landscape at zero temperature for the dipolar 
Heisenberg pyrochlore system and have confirmed the existence of the six-fold degenerate 
Palmer-Chalker ground states. Beginning from one of these classical ground states we 
have employed the techniques of non-interacting spin wave theory. Using bosonic spin- 
deviation operators and an analysis of the quantum mechanical equations of motion we 
have numerically constructed the canonical Bogoliubov transformation and diagonalized 
the spin wave Hamiltonian for the interaction parameter set of Gd 2 Ti 2 07. We have 
obtained the spin wave dispersion spectrum for nearest neighbour exchange and both 
nearest neighbour and long range dipole interactions using the Ewald summation 
technique and confirmed it classically treating each spin as a rotor in its local frame. 
The spin wave excitation spectrum displays four non-degenerate modes and is gapped 
throughout the first zone. Ensuring that a proper bosonic transformation matrix has 
been found for each wavevector k, we have calculated the numerical contribution of 
zero point quantum spin fluctuations to the ground state energy and the reduction of 
the classical staggered moment by performing a discrete average over the first Brillouin 
zone. The partition function was calculated using the fully diagonalized spin wave 
Hamiltonian and we obtained equations for the most relevant thermodynamic relations 
at low temperatures. 

We have confirmed the stability of the PC ground states against quantum 
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fluctuations and it appears that due to the significant anisotropy gap present throughout 
the zone and the large value of the classical spin S = 7/2 of the gadolinium ion, quantum 
mechanical effects are small (~10%) in this system. The fully ordered spin polarization 
is reduced by less than five percent of its classical value, and it should be well described 
by classical methods at low temperature. The physical system retains little 'quantum 
mechanical' knowledge of its fully disordered paramagnetic state after the effects of 
dipole interactions have been included in the Hamiltonian. 

Small perturbations arising from exchange interactions beyond nearest neighbour 
may also be significant in this system, and recent mean field theory studies of Gd 2 Ti 2 07 
with J 2 7^ J3 7^ have identified another possible kagome-like ground state with ordering 
wavevector q = (Ho)- This state has one sublattice with a net moment of zero 
and it appears to be energetically selected over the PC states above some crossing 
temperature T cross [35j. Powder neutron scattering results for Gd 2 Ti 2 7 [12] and recent 
theoretical work by Cepas and Shastry would tend to confirm the existence of this 
1/4-disordered ground state at temperatures below 1 K. In contrast, an experimental 
study of the frustrated pyrochlore material Gd 2 Sn 2 7 shows a single sharp first 
order transition to a long-range ordered magnetic state near 1 K which may be PC in 
character. It is therefore still under debate whether the real low temperature ordered 
phase of these Heisenberg-like dipolar pyrochlore materials is well described by the 
PC manifold. Nevertheless, the stability of the six-fold degenerate classical PC states 
against quantum fluctuations in these systems is of considerable interest. 

Although we have chosen to limit our calculation to one specific material, 
gadolinium titanate, we have presented the calculation in a general form with no real 
dependence on the structure of the lattice, or the spin symmetry. It is therefore hoped 
that the theoretical techniques presented in this study can be employed to provide 
further understanding of the peculiar and interesting properties of highly frustrated 
antiferromagnets on non-Bravais lattices. 
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Appendix A. Numerically Constructed Bogoliubov Transformation 

The equations governing the construction of Q can be derived [HE] by utilizing the vector 
commutator 




(A.l) 
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where X is the column vector of adjoint boson operators and T indicates the usual 
transpose operation. Using this definition we find 

[X, Xt] = g. (A.2) 

where the metric of the transformation is 



g 



1 
-1 



(A.3) 



and 1 is the s x s identity matrix where s is the number of magnetic sublattices. The 
required ortho-normalization condition for the transformation matrix Q is found to be 

QgQ f = g. (A.4) 

The hermitian conjugate of Q can be expressed in terms of its inverse as Qt = gQ _1 g 
and the diagonalization condition is now given by Q^HQ = gQ _1 g = A leading to the 
more familiar eigenproblem 

Q^gHQ = gA. (A.5) 

Defining L = gH and Al = gA we can use a standard numerical package to solve for 
the eigenvectors and eigenvalues of the non-hermitian matrix L, allowing us to write 
QlAlQl 1 = L where Qi_ is the matrix whose columns are the eigenvectors of L subject 
to normalization with respect to the usual identity metric. 

In order to find the proper diagonalizing transformation Q at each fc-vector in the 
first Brillouin zone, three conditions must be met: 

(I) QAlCT 1 = L 

(II) QgQt = g 

(III) Q f HQ = A. 

These can be satisfied simultaneously, provided that we construct a Bogoliubov 
transformation which imposes the full bosonic constraint on our transformed operators. 
This amounts to properly grouping, organizing and then renormalizing the computed 
eigenvectors of L through a block transformation, ensuring that we always remain within 
only one eigenspace at a time [39J. 

Suppose L has m unique eigenvalues, each having degeneracy dj for j = 1, . . . ,m. 
We assume that our eigenvector matrix Ql is organized such that eigenvectors belonging 
to the same eigenspace are grouped together. Our goal is to construct a block diagonal 
transformation matrix dB consisting of m dj x dj blocks which will transform Ql to our 
desired matrix Q via Q = QlcIB. By performing this transformation we have simply 
formed new linear combinations of the eigenvectors within a given eigenspace, and thus 
condition (I) will be satisfied. Condition (II) requires that 

(Q L dB)g(Q L dB) t = g 

and provided that there are only gapped spin wave excitations (no zero eigenvalues), Ql 
is non-singular so we can rearrange this as 

dBgdB f = M, (A.6) 
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where we have defined the block-diagonal matrix M = (Qi_^gQi_) ■ For each block dB 
we must solve an equation of the form 

fcdBjdB} = M ; (A.7) 

where 

fSi = \ , r*. , (A.8) 

I — 1 otherwise 

M J = (Q t Ljg Q Lj )~ 1 (A.9) 

and Qlj is the 2s x dj block of vectors belonging to eigenspace j. An equation of of this 
type can be solved by decomposing Mj = VjDjVj 1 where Vj is the unitary eigenvector 
matrix and Dj is the diagonal eigenvalue matrix corresponding to Mj. Substituting this 
into Equation (jA.7|) we find that 



dB, = VjyfcD-j Vj 1 (A.10) 

is always a solution. Condition (III) is finally satisfied and we have arrived at an 
equation which can be used to numerically construct each block of dB yielding the 
proper transformation matrix Q for each k. 
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